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Abstract 
00 ; 

^^ ■ Employing a rigorous theoretical method for the construction of exact many-electron ground 



"Ti , states we prove that interactions can be employed to tune a bare dispersive band structure such 

that it develops a flat band. Thereby we show that pentagon chain polymers with electron densities 
above half filling may be designed to become ferromagnetic or half metallic. 
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PACS numbers: 71.10.Fd, 71.10.-w, 71.27.+a 



Conducting polymers [l| are a fascinating class of materials with a strikingly wide range 
of applications, e.g., in nanoelectronics |2|, nanooptics |3|, and medicine j^. Many of them 
contain chains of five-membered-rings as a building block. Such pentagon chain polymers 

nn n 

have been explored [5|, |6| and utilized [2| intensively in the past. In particular, polythio- 
phene |7H9| was studied in the search for plastic ferromagnets and, more generally, for 
ferromagnetism in systems made entirely of nonmagnetic elements. The possibility for fer- 
romagnetism in these systems was investigated theoretically 10|, lUi, with a particular focus 
on ferromagnetism due to flat electronic bands arising in odd-membered ring structures 12 |. 
Particular attention was paid to the role of side groups of the pentagon ring, since these may 
cause fiat bands in the band structure. Suwa et al. [13] proposed that ferromagnetism in 
pentagon-chain polymers such as polydimethylaminopyrrole is related to the hybridization 
of narrow a bands with wide vr bands and therefore modelled this polymer by a periodic 
Anderson model. In the latter model the electronic interaction acts site selective within the 
unit cell and the choice of this model [13| was an attempt to account for the different atoms 
on the pentagon chain. 

In this Letter we investigate pentagon chain polymers (see Fig. 1) by a general multi- 
band Hubbard model where the electrons experience local Coulomb interactions on all lattice 
sites. The microscopic parameters are chosen such that they account for the particular 
environment and type of atom in the unit cell of the material; in particular, in our approach 
repulsive on-site interactions are permitted to differ on individual sites. Similarly, we also 
include bond dependent hopping amplitudes. The hopping parameters are not assumed to 
take special values leading to fiat bands in the bare band structure. By contrast we will 
show rigorously that the dispersion of the correlated system may be tuned by the interaction 
to become fiat. Thereby transitions to ferromagnetic states or correlated half- metallic states 
at high electron densities may be induced. We thus prove by exact means the conjecture of 



Brocks et al. IGj that the Coulomb interaction is able to stabilize magnetic order in acene 
and thiophene. 

Our analytic approach proceeds in three steps: the transformation of the Hamiltonian 
into positive semidefinite form, the construction of ground states, and the proof of their 
uniqueness. This technique is independent of the spatial dimension and does not require 
integrability of the model. Previously it was successfully applied to construct exact ground 



states for Hubbard chains with other geometrical structures 17|] and even for the three 




FIG. 1: (Color online) Schematic view of the pentagon-chain polymer polymethylaminotriazole. 
The dotted square indicates the cell presented in detail in Fig. 2. 



dimensional periodic Anderson model 18|. Details of the method are described in Ref. 



17|. 



Since the analytic technique employed here is applicable to a large class of chains we 
discuss it in its general form, and then specify it to a model analysis of the organic pentagon 
chain shown in Fig. 1. The unit cell contains m 



nip + nie > 2 sites where nip is the 



number of sites in the closed polygon as indicated in Fig. 2. In this case mp=5, and rrie = 1 
is the number of sites in the side groups; the number m (here m = 6) also determines 
the number of sublattices. The i-th unit cell then contains the sites i + r„, n < m, with 
ri = and r^+i = a; here |a| is the lattice constant. Altogether the chain consists of A''^ 
cells, where neighboring cells connect through the single point r^+i. Subsequently we use 
periodic boundary conditions and fix the number of electrons to A^ < A^a, where Aa = rnN^ 
is the number of sites. The filling is denoted by p = N/2N\ < 1. In the following ^j, Yl^ 
(°^ X]k' rik i^ momentum representation) mean sums and products, respectively, over the 
Nc cells. 

The Hamiltonian we choose to describe the polymer chain has the form H = Hq + Hjj 
with 

Hq = Y^ Y^ (in,n'4+r„,aCi+r„,,<x + H.C.) + 
(T,i n,n'{n>n') 



a,i n=l 
m 

Hu = YYUnni+r„,tni+r„,i- 

i n=l 
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Here c- ^ creates an electron with spin a at site j, t„_„/ are hopping matrix elements connecting 
the sites i + r,„/ and i + r.„. Furthermore, e„ and [/„ > are on-site potentials and on-site 



Coulomb interactions, respectively, defined at the sites i + r„. The Fourier transform of 
Ci+r„,o- will be denoted by c„,k,o-- We note that the Hamiltonian parameters are arbitrary at 
this point, i.e., they are not chosen to provide fiat bands in the bare band structure. The 
case e„ = for all n leaves the results qualitatively unchanged. 

In the first step, we define m — 1 block operators G'^i^ = Xl^ee- '^a/'^l+ri a with a = 
1, . . . , 771 — 1, i.e., linear superpositions of creation operators acting on blocks B\^a consisting 
of the m sites i + r^ in the unit cell at i; here aa/ are numerical coefficients. Specifically 
for the pentagon cell shown in Fig. 2 we employ 3 three-site blocks (the triangles made of 



sites (1,2,5), (2,3,5), (3,4,5)), and 2 two-site blocks (the site pairs (5,6), (4,7)) [l9|. The 
interaction term Hjj is rewritten in terms of the operators Pn = X]i-^i+r„) where Pj = 
^j.t'^j,-!- ~ ('^j.t + ^j,4.) + 1 is a positive semidefinite operator with eigenvalue zero when there 
is at least one electron on site j. Altogether H — Cg takes the positive semidefinite form 
H — Cg = Hg + Hp^ where 



m— 1 



^^ = E E G.,i,.Gl,^., Hp = J2 UnPn. (2) 



i,(T a=l n=l 

■\m—l -1 Y~^ I |2 



Here Cg = quN - A^cE^=i ^n + 2X1^=1 ^a], Za = Y.t I«a/p5 and qu are constants which 
depend on the parameters entering in H . The transformation of H into the semidefinite 
form shown in Eq. ([2]) requires that the microscopic parameters in Eq. ([T]) fulfill certain 
conditions, i.e., equations connecting the coefficients a^^ and qu to the starting Hamiltonian 
parameters in Eq. ([1]). The solvability of the matching conditions determines the parameter 
space domain T> for which the transformation from Eq. ([1]) to Eq. ([2]) can be performed. 
It is not difficult to show that T> is not strongly restricted by the Hamiltonian parameters 
in Eq. ([T]). In particular, the values of the interaction parameters f/„ can vary over a wide 
range. Details regarding the form of the matching conditions, the domain P, and the block 
operators G]^ j ^ in Eq. (2) are presented in the Appendix. 

Before we construct ground states of Eq. ([2]) above half filling (p > 1/2) we an- 
alyze its effective band structure. This is possible because the anticommutation rela- 
tions for the composite block operators G^ia^ which depend on the interactions f/„, al- 
low us to rewrite the operator Hq as Hg = Hkin + Kg, with a kinetic energy operator 
Hkin = -Y.i,aY.a=i ^l,i,a^o^,h<^ and a constant Kg = ^N^YJ^Zl z^. The kinetic part, 
Hkin, is quadratic in the original fermionic operators c\+rn,a and can hence be diagonalized, 
leading to an effective, interaction dependent band structure. In fact, the dispersion rela- 
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FIG. 2: (Color online) The pentagon cell with nip = 5 and nie = 1- The numbers indicate the 
index n of the site, a is the primitive translation vector, and i + r„ specifies the site position inside 
the cell. 

tions thereby obtained are identical to the energy bands of Hq in Eq. ( ITb|) . but the on-site 
potentials e„ are now replaced by the renormalized energies 

ef = en + f/n - qu- (3) 

It is this renormalization which can lead to an effective upper fiat band. In the ground state 
|\I'g) of H with {Hq + Hp)\'^g) = (see below), this flatness is unaffected by the presence 
of Hp in Eq. ^ since Hp\^g) = 0. 

Thus we find the very remarkable result that a dispersive band structure of noninteracting 
electrons can be tuned by an interaction to yield an effective upper fiat band of the inter- 
acting, many-electron system; an example is given Fig. [3] for a selected parameter set. The 
effectively flat band is half filled for the total number of particles N = 2Nx — Nc = N*, and 
is more than half filled for A^ > A^*. Such an interaction induced upper fiat band is possible 
only if the local interactions f/„ differ on at least one site in the unit cell. We note that 
properties of the exact ground state can only be deduced for the states in the upper band, 
the physics of the lower bands is not accessible by the here applied method (as indicated by 
the question mark in Fig. 3b). 

The ground state for N = N*: In this case the ground state of Eq. (|2]) has the form 

(T 

where |0) is the vacuum state, Gj. = 11111^=1 ^mio-' ^^^ ^^^ operator F'l" = Hi^l+r o- 
introduces one electron with spin a in each unit cell. Since Gj. creates {m — l)Nc electrons 



with spin a, the state |\I/g(A^*)) contains N^r = mNc electrons with spin a. Therefore there 
is one a electron on each site. Consequently all a electrons are localized, and only the —a 
electrons are mobile. Therefore Eq. (^ may be identically rewritten as 

m m—1 

\^,{N*)) = n[(n^u,j(n ^u-mo)- (s) 

i n=l a=l 

This state describes a half-metal, i.e., a non-saturated ferromagnet with total spin S = Nc/2. 
Eq. (4) is indeed the ground state since Hd'^'giN*)) = and Hpl'^g^N*)) = 0, where the 
former relation is due (Gaia)'^ — 0' while the latter is a consequence of N^ = Na- Since 
every site is occupied by an electron with spin a that part of the wave function which 
describes the —a electrons is equivalent to a Slater determinant. 

The spatial extension of the electrons with spin —a is obtained from the long-distance 
(r — 7- oo) behavior of the ground-state expectation value of the hopping term ri(r) = 
{^ giN*)\{cl^j. _^Ci^r„+r-a + H.c.)\'^ g{N*)) for arbitrary n = 1, . . . , m. Explicit calculations 
yield an exponential decay of ri(r) in the thermodynamic limit {Nc — > oo). Hence the 
ferromagnetic state Eq. (5) is localized. Apart from the trivial (25* + 1) degeneracy related 



to the orientation of the total spin, where S = Sf'-'"'^ = Nc/2, the ground state [20| is unique; 
this was proved by us using the technique presented in detail in Ref. |l7l |. 

The ground state for N > N*: Here we restrict ourselves to the 5*2 = 3^'^°'^ sector. If we 
add N electrons to the system these electrons can occupy only —a spin states. The ground 
state then has the form 

\^g{N* + N))=Q^^\^giN*)), (6) 

where the operator QJy = [117=1 ^I k -a] ^^ ^ product of N arbitrary, but different, c|jj^_^ 
operators. That Eq. ([6]) is indeed a ground state follows from the fact that G]^ ; ^ anticom- 
mutes with the fermionic creation operators cl^ ^ _^ and hence also with the operators in QU. 
Due to the free electrons with spin —a introduced by Qjy the ground state Eq. ([6]) contains 
also electrons in plane- wave like states. Furthermore, since Eg = Cg, one finds a vanishing 
charge excitation gap 6fi = Eg{N + 1) - 2Eg{N) + Eg{N - 1) = ioi N > 1. 

In order to verify the extended character of the ground state Eq. ([6]) we calculate again the 
expectation value of the hopping term ri(r) = {'^g{N* + N)\cl^^^_^Ci+r2+r-a + H.c.\'^g{N* + 
N)) specifically for N = 1 and Qi = cL (i) _ , where k'-^^ is the arbitrary momentum of the 
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FIG. 3: (Color online) Dispersion of the pentagon polymer chain, a) Bare band structure (for 
details see the Appendix) for hopping parameters tc = t^j = 0.5, th = ^3,2 = —1-1, tf = t^^ = 1.2, 
and on-site potentials ei = €4 = —2.5, €2 = £3 = —2.0, £5 = —2.1, eg = —2.1, using the site notation 
from Fig. 2. b) In the interacting case with the local interactions Ui = U4 = 0.18, U2 = U3 = 0.35, 
U5 = Uq = 0.03 and the same Hq parameters as in a). The upper band, indicated by a thick line, 
is determined by Eq. ([3]). Since E{k) is even in A: = k • a, only k E [0, vr] is shown. The question 
mark in the shaded area indicates that the exact behavior in this region is not known. All energies 
are in units of t = ^2,1 = ii,5 = ^5,4 = ^,3, with t > 0. 

electron added above N = N*. In the thermodynamic limit of the pentagon chain the result 
is 

ri(r) = ro(r)(l-A(k«)/fi(k«)). 



(7) 



Here ro(r) is the plane- wave result for a free electron in a Bloch state, A(k) = ^41+^2 cos k-a, 
and B(k) = Bi + i?2Cosk ■ a > holds for all k with constants Ai, A2 and Bi, B2 (for 
details see the Appendix). This shows that the localization length is indeed infinite. 

For pentagon chains without external links (i.e., when the sites i + rg and i + r7 in Fig. 
2 are missing) the same solutions are found, but the regions of parameter space where they 
exist are shifted. 

We emphasize that in the case of Mielke-Tasaki ferromagnetism, i.e., flat-band ferro- 
magnetism in a half-filled lowest fiat band 12|, both the fiat band and the connectivity 
conditions (the overlap of the local Wannier functions) necessary for the emergence of this 
type of ferromagnetism, result from the bare band structure determined by Hq. By con- 



trast, we have shown here that the interaction may be employed to tune a fully dispersive 
bare band structure to become partially flat. The results apply to the high-density region 
where Brocks et al. conjectured the presence of strong correlations in acene and thiophene 
organic molecular crystals and the stabilization of magnetic phases |il&]. Using the Hq pa- 
rameters from Fig. 3, the matching conditions provide for the upper (bare) band a width 
of W/t = 0.15, with 0.2 < Un/W < 2.3. For other parameters in Hq the parameter space 
domain V (see the Appendix) allows for even higher values of Un/W. 

Regarding the experimental realization of the ferromagnetic and the half-metallic states 

derived above we note that the required electron doping of the pentagon chains can, in 

principle, be achieved |2l| by two means: first by raising the Fermi level by selecting appro- 

)riate side groups, and second by field-effect doping in a double-layer transistor structure 

22 1 ■ Specifically for polythiophene with an estimated density of 10^^ pentagon rings/cm^ 

231], electron densities of the order 10^^ — 10^^ carriers/cm^ are required. As already verified 

these electron densities are experimentally achievable |24l |. 

In summary, by employing a rigorous analytic method we have constructed exact ground 
states for a multiorbital pentagon Hubbard chain. The ferromagnetism and the half- 
metallicity of the derived solutions originate from an unexpected mechanism in multiorbital 
polygon chains with different site-dependent Coulomb interaction strengths. For high elec- 
tron densities with the top band at least half-filled the interactions are capable of turning 
this dispersive band into an effectively fiat band in an extended parameter region. The ob- 
tained solutions therfore point to a new route for the design of ferromagnetic pentagon-chain 
polymers. 
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APPENDIX 

1. DETAILS REGARDING THE TRANSFORMED HAMILTONIAN, EQ. (2), FOR 
THE PENTAGON CHAIN 

a) The parameters of the transformed Hamiltonian, Eq.(2), obey the matching conditions 

m— 1 

Y,^^ = 4Qi + Ql/\th\ + 2tyQ, + 2(|U| + |t,|) + QI + t}/Ql, 

a=l 

qu = (l/2)[(f/i + ei + |t,|) + (f/2 + 62 + \th\) + {[(f/2 + 62 + \th\) - (f/i + ei + M? + ^t'Y% 

Qi = Qu - U2 - £2 - \th\, Q2 = qu - Ui - ei - |tc|, 

Qs = \tf\VK\/[\th\{qu - f/5 - es) - {qu - U2 - £2)' + tlf', (Al) 

where Qi,Q2,Q3 > 0, and the hopping amphtudes t,th,tc,tf are defined in the caption of 
Fig. 3. 

b) The parameter space domain V for which the transformation from Eq.(l) to Eq.(2) is 
vahd is determined by the relations: 

t,, < 0, Z = {qu - Ql) > 66, 

W = qu- [{qu - f/2 - 62)' - tl]/\t,\ > 65, 

W-e,>U,>0, U^ = Z-e^. {A2) 

10 



c) The block operators entering in the transformed Hamiltonian, Eq.(2), are given by the 
relations 

G^2,i,o- = —tsbi,2,3 + QlCi+r5,a/ts, 

G3,i,a = (lsbi,5,3 — t/qsCi+r4,(T, 

Cf4,i,o- = QsCi+rs.o- — ^//v3Ci+r5,o-, 

G 5,i,a = \tc\ [Ci+r4,(T ~ ^cCi+a,cr/Kc|J) (^"j) 

where the notation qs = Q^"^ , ts = |t/l|^/^ &i,„,m = Zlp=„,mCi+rp,a was used. 

2. DISPERSION OF THE BARE BANDS OF THE PENTAGON CHAIN PRE- 
SENTED IN FIG. 3A 

The energy of the bare bands presented in Fig. 3a, e = E^{k), u < m, is provided by the 
equation 

2t,t'[{e, - e)T, - t^Tj] cosk + Tf{tltl - {e, - eftl + [{e, - e)(e2 - e) - t^f - (ei - e)Hl} 
+2(e6 - e)e{t\e2 + t^ - e) - (ei - e)Tn} = 0, (A4) 

where Th = (e2 - e)^ - t|, Tj = (ee - e)(e5 - e) - tj, and A; = k ■ a. 

3. DETAILS REGARDING THE EXPECTATION VALUE OF THE HOPPING 
TERM ri(r) 

The parameters entering in the ri(r) expression from Eq.(7) are given by the relations 
i2ro,n2/oij- I I /n ^ I n-^2fr^ i i-* l^ i /n /^2 , o/n2N , /rM , /n2j.2\ 



A, = 2\t,\{Qi[2Qii2m + go + 2t\Q, + |t„|) + Q,it' + 2Qi) + (Qt + Qit')/^] 
4/ 



+tyQl[t'{\tr.\ + go + Ql - {Qt + Qlt')/K\]}, 



A2 = 2t,tWtr,\t})/Ql - 2Ql{Q, + |t. 



B2 = 2t\[{\tn\t})/Ql - Qli\tn\ + 2Q,)], (A5) 



B^ = 2|t,|{g^[(2|t„| + Qi/K\){Qi + n + it' + 2Qi){\tn\ + 2Q,)] + {K\tHj)/Q',} 

and one has Bi + B2 cos{k) > for all A; = k ■ a. 

11 



